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SUMMARY 

A method recently developed for determining the steady flow of a 
nonviscous compressible fluid along a relative stream surface between 
adjacent blades in a turbomachine was applied to investigate the sub- 
sonic through flow in a single-stage axial -flow gas turbine. A free- 
vortex type of variation in tangential velocity was prescribed along the 
stream surface. Cylindrical bounding walls were specified in order to 
avoid radial flow at the walls. The flow variations on the stream sur- 
face for incompressible and compressible flows were obtained by using 
the relaxation method with hand computation and the matrix method on both 
an IBM Card Programmed Electronic Calculator and a UNTVAC . 

In all solutions considered in this investigation, convergence was 
obtained without difficulty. A comparison between the relaxation and 
the matrix methods showed that more accurate results were obtained with 
the matrix method in a shorter interval of time. The results of these 
accurate calculations provide a basis for evaluation of simpler, more 
approximate methods for computing subsonic through flow in turbines. 

Considerable radial flow was obtained for both incompressible and 
compressible flows because of the radial twist of the stream surface 
required by the prescribed velocity diagram and the compressibility of 
the gas ; the radial twist of the stream surface and the compressibility 
of the gas had equally important effects, and the nonlinear nature of 
the equations defining the flow was quite evident . The shape of the 
stream surface was found to be sensitive to the axial position of the 
radial element of the stream surface in the stator. When the radial ele- 
ment of the stream surface in the rotor was near the midaxial position in 
the rotor, a large negative gradient of axial velocity was observed in 
all cases ahead of the rotor. 
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INTRODUCTION 

In the design of gas turbines having relatively long blades in the 
radial direction, a first approximation of the radial variation of the 
state of the gas is commonly obtained by assuming the gas to be nonvis- 
cous and to flow on cylindrical surfaces. Axial symmetry is assumed in 
such one -dimensional- solutions for the radial variation in the gas state 
upstream and downstream of a blade row and they are usually referred to 
as the "simple-radial-equilibrium" or "simplified-radial-equilibrium" 
solutions (see references 1 and 2 for examples). Even turbines designed 
from a free-vortex velocity diagram, however, have considerable radial 
displacement of gas particles across the blade row due to the compressi- 
bility of the gas; this displacement, of course, violates the assumption 
of flow on cylindrical surfaces. The radial displacement of gas particles 
across the blade row and the curvature in the streamline caused by the 
radial motion have a significant effect on the radial variation in the 
gas state (see reference 2). 

Several through-flow methods have been developed to consider both 
radial and axial variations in flow conditions (see, for example, refer- 
ences 2 and 3) . In these methods an infinite number of blades and a 
nonviscous fluid are assumed. The interpretation of such a solution, 
its modification for any finite number of thick blades, and it's exten- 
sion to a complete three-dimensional solution for nonviscous fluids are 
given in reference 4. In actual turbines, the motion of the gas is 
further complicated by the secondary flow caused by the boundary layers 
along the hub and casing walls (for example, reference 5). Nevertheless, 
it is believed that the detailed flow analysis based on an assumption of 
a nonviscous gas will permit a clearer understanding of the individual 
contributions to the complete flow and help to explain the development 
of the viscous boundary layer along the walls. 

The theoretical method of reference 4 for obtaining a complete 
three-dimensional solution is based on an appropriate combination of a 
number of mathematically two-dimensional flows on two kinds of relative 
stream surface. The first kind of relative stream surface extends from 
the suction surface of one blade to the pressure surface of the adjacent 
blade and deviates, in general, from a surface of revolution about the 
axis of the turbomachine; in reference 4 this surface is called Sp. The 

second kind of relative stream surface extends from hub to casing and 
roughly approximates the shape of the mean camber surface of the blades; 
in reference 4 this surface is called S 2 . The equations defining the 

flow on these two kinds of stream surface sire, however, similar and the 
methods of successive approximation used in their solution sire essen- 
tially the same. The mechanics of obtaining a numerical solution of 
either of these two-dimensional problems by the methods of reference 4 
has not previously been investigated. If convergence of one of these 
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two-dimensional solutions can be obtained without difficulty, the other 
solution should then converge as well. The method of reference 4 will 
then very likely be practical for analyzing three-dimensional flows of 
compressible nonviscous fluids in turbomachines having finite numbers 
of thick blades and arbitrary hub and casing shapes. 

The primary objective of this investigation made at the NACA Lewis 
laboratory is to study in detail the technique of solving for the sub- 
sonic flow along a relative stream surface of the second kind (S2) by 

using the methods of references 4 and 6. At the same time, of course, 
the effects of compressibility and radial twist of the stream surface 
are obtained; the examination of these effects therefore constitutes 
the secondary purpose of this investigation. If highly accurate solu- 
tions can be obtained in this way, these solutions will provide a basis 
for evaluation of other simpler but more approximate solutions for sub- 
sonic through flow in turbines. 

For these calculations, a single-stage axial-flow turbine was 
selected as the machine type to be analyzed. The annulus walls were 
chosen to' be cylindrical and to have a hub-tip radius ratio of 0.6. 

The radial distribution of the tangential component of velocity was of 

the free-vortex type. The turbine-work parameter -AH/U^-^ was 0.96 
where -AH represents the stagnation-enthalpy drop and U^, the blade - 

tip speed. The tangential component of velocity at the rotor exit was 
assigned equal to zero. For a compressible fluid, the assigned condi- 
tions resulted in an absolute Mach number of 0.99 at the root and exit 
of the stator and a relative Mach number of about 0.70 at the root and 
entrance of the rotor. The effects of radial blade force were varied 
by changing (l) the axial distribution of aerodynamic loading and 
(2) the axial location of the radial element of the stream surfaces in 
the stator. 

The following three methods were employed in making the calcula- 
tions : 

(1) Relaxation method with the use of a hand-operated desk calcu- 

lator 

(2) Matrix method on an IBM Card Programmed Electronic Calculator 

(3) Matrix method on a UNIVAC 

The matrix factors obtained in the matrix solutions can be used for 
similar calculations for compressors and turbines which have a constant 
hub -tip radius ratio of 0.6. 
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SYMBOLS 

The following symbols are used in this report. (Bold symbols indi 
cate vectors on the stream surface; nonbold symbols indicate the compo- 
nents of these vectors.) 

a velocity of sound 

B variable defined by equation (2c) 

V differentiation coefficients in equation (26) used to 

4 3 multiply function value at grid point j to give m 

derivative at grid point i using polynomial of 
fourth degree 

[c],Ce], 

Col, Cl 3, 

_D * 

Dt 

F 

H 


h 


L square matrices 

hi 

differentiation with respect to time following motion of 
gas particle 

vector having dimensions of force per unit mass of gas; 
defined by equation (6) 

total, or stagnation, enthalpy per unit mass of gas, 
h + H* 4 

enthalpy per unit mass of gas 


i,j,k grid points 

M Mach number 

N nonhomogeneous term of principal equation as given in 

equation (5) 

n unit vector normal to stream surface 

p static pressure of gas 

any dependent variable 
R gas constant 

r radial distance , r*r^ 
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s 

entropy per unit mass of gas, s*R 

T 

temperature of gas 

V 

t 

time 

U 

velocity of blade at radius r 

5f v 

absolute gas velocity, V*Ut 

X W 

relative gas velocity, w*U-fc 

Z 

axial distance, z*r t 


column matrices 

r 

ratio of specific heats, 1.33 

5 r* 5 z 

grid spacing in r- and z-directions, respectively; 
^r r t > ^z r t 

P 

mass density, p*P T i 

E,® 

general variables used in density table 

<P 

angular position 

f 

* *? 

stream function, \|r PT,iU-t r t 

CD 

angular speed of blade 

Subscripts : 


e 

exit 

h 

hub 

i 

inlet 

A 

grid point 

o 

refers to position where stream surface has a radial ele- 
ment or where F r = 0 

r 

radial component 

T 

total, or stagnation, state 
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t tip 

u tangential, or circumferential, component 

z axial component 

Superscripts : 
i,j,k grid points 
* dimensionless value 


EQUATIONS GOVERNING FLOW ON A RELATIVE STREAM SURFACE 
BETWEEN TWO ADJACENT BLADES 

The present report is concerned with the inverse solution (design 
problem) for the steady compressible flow on a relative stream surface 
about midway between two adjacent blades (see fig. l) . The shape of 
this stream surface is not known in advance but the variation in tangen- 
tial velocity of the gas and the position of the radial element of the 
surface are prescribed in the design. The shape of the stream surface 
as well as the state of the gas flowing along the surface is described 
by the two independent variables r and z. 

In reference 4, the following continuity equation for steady flow 
on the stream surface is obtained: 


8(rBpW r ) 3(rBpW z ) 

3r 3z 


( 1 ) 


In equation (l) the bold partial derivative sign denotes the rate of 
change of the dependent variable on the stream surface with respect to 
the independent variable and is related to the ordinary partial deriva- 
tive with respect to the coordinates r, cp, and z as follows: 

2a = |a . H i |a ( 2a ) 

3r dr n^ r ocp 

§£ = ^ i (2b) 

9 z oz n^ r ocp 


The angular variation of the variable is thus implicitly included 
although its value can be calculated only after the shape of the stream 
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surface is obtained in the solution. The variable B in equation (l) 
is related to the angular variation of the gas velocities and the shape 
of the stream surface by 


B f /*r dW r i n z dW z \ 

B \^W + T^ + ^W/ dt 


(2c) 


This variable B can also be interpreted as a variable angular thick- 
ness of a thin stream sheet whose mean surface is the stream surface 
shown on figure 1 (reference 4) . Because B along the mean stream sur- 
face is very closely related to the ratio of circumferential channel 
width to pitch (reference 7), it can be approximately estimated in the 
design by considering desirable axial and radial variations in blade 
thickness required to provide, for example, adequate strength and a 
cooling passage of sufficient size; with these variations taken into 
account, B can then be used in the design calculation of the mean 
stream surface. For the limiting case of an infinite number of blades 
of zero thickness, B becomes a constant and is taken equal to 1. 

Equation (l) is the necessary and sufficient condition that a 
stream function \|r exists with 


If - rB F*z 

(3) 

= -rBpW„ 

az 

(4) 


With the tangential velocity of the gas specified* the subsonic 
flow of the gas is obtained by the solution of i|r in the following 
principal equation, which is obtained from the equation of motion in 
the radial direction with the use of relations (3) and (4) : 

£ ft l ai|f a 2 t _ TJ 

3r 2 r 9r 3Z 2 

where 


1 

a (Bp) 

a* 

a(B P ) 

3\jr 

_ (rBp) 2 


3(V) 

_ 3H _j_ ip as p 

Bp 

ar 

ar 

az 

az_ 

a^/ar 

_ r 

3r 

ar 3r r _ 


( 5 ) 
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In equation (5), F r is the radial component of a vector F, which is 
defined as follows : 

F . - - — n (6) 

n u r\^3cp ocp J n u rp dcp 

For the limiting case of an infinite number of blades, F becomes the 
blade force. For the evaluation of F r , which represents the influence 

of the radial twist of the stream surface and of the circumferential 
pressure gradient on the gas flow, the other components of F are 
first computed by the equation of motion in the tangential and axial 
directions : 



at a(y) + 9 jr a(y) 

dz ar ar az _ 


(?) 



The component F r is then obtained by the use of the following equa- 
tion, which is derived from the integrability condition, which insures 
that the stream surface to be obtained is a continuous integral sur- 
face : 


r> z 



The radial derivative of entropy s in equation (5) is determined 
from its radial distribution at the inlet and the following condition 
of the constancy of s along a streamline on the stream surface for 
reversible adiabatic flow: 


= 0 
Dt 


( 10 ) 
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The radial derivative of H in equation (5) is determined from the 
inlet variation and the following equation: 


DH D (V) 
Dt ~ “ Dt 


(ID 


Equation (ll) is obtained (reference 4) by the use of the equation of 
motion on the stream surface and the following equation expressing the 
orthogonal relation between the resultant relative velocity and the 
surface normal n or its parallel vector F : 

WrF r + W U F U + W Z F Z =0 (12) 


Because the three equations of motion have already been employed in the 
solution, there is only one more independent relation in equations (ll) 
and (12). In the following, equation (ll) is used and is considered to 
represent the orthogonal relation (12). 


The variation of density included in N in equation (5) is deter- 
mined from the ^-derivatives by use of the following relation between 
density, enthalpy, and entropy, between the local static condition and 
the total condition at the inlet : 



PRESCRIBED DESIGN CONDITIONS 

In the present study of the through flow in a gas turbine, the 
effects of some of the design variables are considered. Cylindrical 
bounding walls are specified in order to avoid radial flow at the 
walls. The meridional section of the turbine is shown in figure 2. 

A hub-tip radius ratio of 0.6 and a blade aspect ratio of 2.67 (which 
is based on axial chord and corresponds to the blade-row aspect ratio 
of 2 used in reference 2) are chosen in order to compare some of the 
results with those previously obtained in reference 2 by an approxi- 
mate method. 
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Flow conditions were computed for six sets of assigned conditions. 
These sets of conditions are designated cases A, B, C, D, E, and F 
and are summarized in the following table : 


Case 

Fluid 

Loading 

distri- 

bution 

Axial location 
in stator of 
radial element 
of stream 
surface 

* 

z o 

A 

Incompressible 

Uniform 

0.0675 

B 

Incompr e s s ib le 

Nonuniform 

.0625 

C 

Compressible 

Nonuniform 

(F r neglected) 

D 

Compressible 

Nonuniform 

.0625 

E 

Compressible 

Nonuniform 

0 

F 

Compressible 

Nonuniform 

.1500 


The prescribed variation of tangential velocity, or the angular 
momentum per unit mass of gas V u r*, is such that at a constant-z plane, 
V*r* on the stream surface is constant with respect to r *; and at all 
fixed values of r* , the variation of V u r on the stream surface with 
respect to z* is as shown in figure 3. Two kinds of variation with 
respect to z* are considered. The dashed line shows a linear varia- 
tion, and the solid line shows a composite variation in which a constant 
rate is maintained for the first half of the blade chord and a rate 
linearly decreasing to zero is used for the second half. These two var- 
iations are called uniform and nonuniform loading, respectively. In 
both cases the total change of V u r divided by U^r^ (or -AH/U^ 2 ) 

across the blade row is 0.96, which is used in reference 2. The expres- 
sions for the dimensionless specific angular momentum are as follows: 

(a) Uniform loading 

Stator: 0 < z * < 0.15, 


a(V* r *) _ 0.96 
9 z * °- 15 


V = 


0.96 


(14) 


0.15 


(15) 


2402 




NACA TN 2750 


11 


Rotor: 0.20 < z* < 0.35, 


a K r *) _ o 


az 


96 
0.15 


<r* = ^|| (0.35 


0.15 

(t>) Nonunlform loading 

Stator: 0< z* < 0.075, 

r * * 


**) 


(16) 

(17) 


Rotor : 


9 (V ) 4 0.96 


az 


3 0.15 


* * _ 4 0-96 * 

u r 3 0.15 


(18) 

(19) 


0.075 < z* <0.15, 


9 (V* r *) 8 0.96 / z* \ 

az * 3 0.15 \ " 0.15/ 


( 20 ) 




- 0.32 + ± X 0.96 


— ~(— ) 2 ] 

0.15 \0 . 15 / 


( 21 ) 


0.20 < z* < 0.275, 


8 ( v u r *) _ 4 0.96 


az 


3 0.15 


V*r* = 0.96 - £ (z* - 0.2) 


3 0.15 


( 22 ) 


(23) 


0.275 < z* < 0.35, 


3(vM) 


az 


8 0 
3 0 


^■h - 

.15 \ x 


z - 0.2 

0.15 


( 24 ) 
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vjr* = 1.28 - £ X 0.96 

For a given axial distribution of tangential velocity, a change in 
the axial location of the radial element of the stream surface will 
alter the radial and axial distributions of both axial and radial velo- 
cities and thereby change the shape of the stream surface. In the 
rotor, the axial location of the radial element of the relative stream 
surface is 45 percent of the axial blade chord from the leading edge 

(z* = 0.2675) for case A and 41.7 percent (z* = 0.2625) for cases B, 

D, E, and F. 

With this prescribed axial variation in V u r, the radial deriva- 
tive of V u r contained in the equations given in the preceding sec- 
tion drops out. Also, for the present investigation, the inlet flow 
is considered to be uniform in entropy and total enthalpy. Then, for 
adiabatic frictionless flow, the radial and axial derivatives of s 
vanish. The radial derivative of H also vanishes for the specified 
inlet condition and V u r . 

For the present investigation, B is taken to be a constant (a 
value of 1 is used in the calculation) . In a sense , the solutions 
thus obtained do not depend on any particular blade configurations. 

But the solution is correct only for those bladings whose geometrical 
configuration is such that the angular thickness of the mean stream 
sheet is essentially constant. It also gives the limiting solution 
for an infinite number of infinitely thin blades. To interpret the 
results obtained as this limiting solution, the prescribed variation 
of V u r corresponds exactly to the free-vortex type with S(V u r)/Sr 

equal to zero. For the general interpretation of the results obtained 
as the solution for the flow along a mean stream surface (subject to 
the assumed constant value of B), the prescribed condition that 
a(Y u r)/ar be equal to zero does not give d(v u r)/c)r equal to zero. 

The inlet flow of the turbine example given in reference 2 is 
used in the present investigation and is as follows : 


v 



M. = 0.308 

l 



✓ 
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For this Mach number. 


p* = — — = 0.95033 
1 0T,i 


Hi* 



12.546 


METHOD OF SOLUTION 

The principal equation to be solved Is equation (5). This partial 
differential equation is first replaced by a number of finite -difference 
equations representing the differential equation at a number of grid 
points covering the domain. Because of the nonlinear nature of the 
problem, these equations are solved by the general method of successive 
approximations. In each cycle of calculation, the nonhomogeneous 
term N is evaluated by employing any appropriate approximate solution 
at the start and by using equations (7) to (ll) and equation (13) 
together with the flow variation obtained in the preceding cycle; the 
results are then taken as given values in the solution of >}r in the 
following cycle. The solution of \|r from the finite -difference form 
of the principal equation is obtained by the relaxation method (refer- 
ence 8) in the modified form as given in reference 6 and by the matrix 
method discussed in reference 6. The accuracy of the solution depends 
on the accuracy of the finite -difference representation of the partial 
differential equation, the size of the residual left in the solution 
of \|r, and the number of cycles completed for convergence. 


Choice of Grid System and Degree of Polynomial Representation 

In these calculations, a single grid size was used for both the 
relaxation and the matrix calculations. The results of reference 2 
were useful in selection of the grid size. From the results obtained 
in reference 2, the stream function \|r is expected to increase 
smoothly with respect to radius at constant-z planes and to vary 
approximately as a sine curve with respect to z at constant r-values. 
With the necessity of covering a large domain in order to satisfy the 
boundary conditions given far upstream of the stator and far downstream 
of the rotor, and with such a smooth variation of \|r over the domain, 
the use of a fourth-degree polynomial representation rather than of the 
usual second-degree one is suggested so that the number of grid points 
may be reduced. From experience gained while obtaining relaxation 
solutions, a final grid size of 0.05 in the radial direction and 
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0.025 rt in the axial direction in the meridional plane are used 
(fig. 4). These grid points on the meridional plane are for reference, 
or recording, purposes only. The variables involved in the solution 
are those on the stream surface (fig. 4). 

The use of this grid size gives seven radial and axial stations 
on the stream surface between the hub and the casing and across each 
blade row. With the present variation of \|r in the radial direction, 
sufficiently accurate results are expected for radial derivatives. The 
accuracy of the axial derivatives is analyzed for a simple sinusoidal 
variation in the axial direction with a period of 20 5 Z , and the fourth- 
degree polynomial representation is found to give first- and second- 
order derivatives accurate to within 0.02 percent. 

For simplicity (in order to give a uniform formula in the 
z-direction) , the same grid size in the z-direction is used for the 
entire domain. It is found in the relaxation solution that the radial 
distribution of i|r had no axial variation in the first five signifi- 
cant figures after eight of these z-stations either upstream of the 
stator or downstream of the rotor. The matrix factors for matrix solu- 
tions therefore cover a range of z* varying from -0.5 to 0.85, which 
includes 10 stations each way upstream of the stator and downstream of 
the rotor. The order of matrices is thus 7x55 = 385, which is also 
the total number of interior grid points. 

At the first and the last z-stations and a few stations nearby, 
sufficiently accurate z-derivatives can be obtained by the use of a 
three -point differentiation formula. For simplicity in setting up the 
matrix factorization, however, the same central -point fourth-degree 
differentiation formulas are used for the entire domain. The use of 
these formulas means that the same boundary values of \|r are used for 
the two stations outside the first (z* = - 0.5) and the last (z* = 0.85) 
z-stations. 


Finite -Difference Form of Principal Equation 

With the grid sizes of 0.05 and 0.025 chosen for 6* and 6*, 

respectively, the differentiation coefficients for the first- and the 
second-order derivatives are computed. If these coefficients at grid 

point i are denoted by ^ and respectively, the finite - 

difference form of the principal equation at any grid point i becomes 
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(26) 


where \jr J and i|r k denote the values of t on the stream surface 


corresponding to the grid points along constant-z and constant-r lines 
on the meridional plane, respectively (fig. 4). 


The value of i|r along the hub is chosen to be zero. The value of 
a dimensionless \|r along the casing is chosen as follows: At sta- 

tion i-i, with the use of equation (3), 


which is a constant along the intersecting curve of the stream surface 
and the casing. 

At the inlet station i-i, which is one 8 z -distance upstream of 
the first station (z* = -0.5), and the exit station e-e, which is one 
5 Z -distance downstream of the last station (z* = 0.85) (see fig. 2), 
the radial distribution of i|r* is confuted as follows: 


Calculation of Boundary Values of \(r 




For equal to 1 and the chosen values of p* and V* 


= 0.19767 


(r) = i|r* (r) = 0.19767 



(27) 
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This simple variation results from the fact that with no "tangential and 
radial velocities of the gas and with uniform values in H and s, the 
axial velocity and density are also uniform at "both these stations, 
qq^is radial variation in i|r is maintained constant from cycle to cycle 
and the problem is therefore treated as a boundary-value problem of the 
first kind. 

If the nonhomogeneous term N is everywhere equal to zero, the 
differential equation (5) would be satisfied by the ifr -function as given 
by equation (27), which means that the gas flows on cylindrical sur- 
faces. The present problem is then essentially to determine the change 
of i(r from this simple distribution due to a certain distribution of 
N in the domain resulting from the nonzero values of density deriva- 
tives and F r . 


Calculation of Nonhomogeneous Term 

Calculation of density by the use of general table . - The defini- 
tion of N in the principal equation (5) shows N to consist of two 
terms for the present investigation. One is connected with the com- 
pressibility of the gas, and the other, with the tangential pressure 
gradient of the gas and the radial twist of the stream surface. The 
first term in N vanishes for incompressible flow, and its evaluation 
for compressible flow is as follows: 

In each cycle of computation starting with a given variation of 
\|r, its derivatives with respect to r and z are first obtained by 
numerical differentiation. These derivatives are squared, added, and 
used in the following formula to obtain a function $ 


_ t+1 



The last two factors in the right side of the preceding equation are 
obtained from the given inlet condition and equation (ll). From the 
value of , a value of L is read from table l(b) given in refer- 
ence 4 (with i = 4/3) , and p* is obtained by 


P* = i 57 




Hi 


2 


(29) 
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After p* is evaluated, its derivatives are computed and combined 
with the i|r- derivatives to form the first term in N. 

Calculation of F r by equation (9) and approximate formulas. - 

The second term in N involves the use of p and 3\|r/9r based on 
the new values of i|r in each cycle and the evaluation of F r in the 
cycle. In order to evaluate Fr, equations (7) to (9) are used. The 
computation of F u r by equation (7) is relatively simple (the second 
term drops out for the present problem) . For the computation of F z 

of equation (8), the last factor of the last term 


~§!i _ 

iat 9 2 t 

_ 1 / 9p 9jr + 9p 

_9r 2 

r ar + 2 " 

3z 

p V ar ar az az ) 


is available in the solution by the relaxation method because it is 
involved in the computation of the residual. In the matrix solution, 
however, this factor is not available. In such a case, the last term 
is replaced by the following expression through the use of equation (l) : 


F 
r r 

in which the value of F r of the previous cycle is used. After F z 
is computed, it is divided by F u r and differentiated with respect to 
z. With z Q chosen between two grid points, the formula given in 
reference 9 is used for the integration to determine F r . 

The following approximate expression for F r is obtained from 
equations (7) to (9) by assuming af/ar constant and neglecting small 
terms : 


9i|r/az 

aV/ar 


F„ = .£. 


2 8(y) 

3 az 


Pz 


(V u r) dz 


(30) 


When equations (14) to (25) are substituted into the preceding 
equation, the following expressions for dimensionless F r are obtained: 
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(a) Uniform loading 


Stator : 


Rotor : 



(31) 


z 


* 2 -0.2 2 j 


(32) 


(b) Ronuniform loading 

Stator : 0<z*<0.075, 


2 4 0.96 

r *3 3 0.15 


„ *2 *2 
4 0.96 z z o 

3 0.15 2 . 


(33) 


0.075 < z* < 0.15, 


F 


* 

r 


2 8 0.96 / z *\ 

r *3 3 0.15 \ " 0.15/ 


4 0.96 (° •' 0752 "Zo 2 
3 0.15 \ 2 


0.32 (z*-0.075) +| (z* 2 -0.075 2 ) - 


1 z* 3 -0.075 3 
9 0.15 


(34) 


Rotor: 0.2 < z* <0.275, 



2 4 0.96 

*3 3 0.15 

r ^ 


0.96(z*-z*) - 


4 0.96 ( z *2 ) 2 -( z J-0 .2 ) 2 
3 0.15 ~ 


(35) 


0.275 < z* <0.35, 
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Fr 


2 8 0.96 L _ z*-0.2\ 

r *3 3 0.15 V 0.15 ) 


0.96(0.275- z$) 


4 0.96 0 . 075 2 -( zJ - 0 . 2) 2 
3 0.15 2 


+ 1.28 (z*-z*) - 



Solution of Principal Equation by Relaxation Method 

Computation of the residual . - After the nonhomogeneous term 
is obtained at each grid point, it is subtracted from the left side of 
equation (26). The difference obtained is the residual at the grid 
point i. 

Reduction of residuals . - The coefficients used to relax residuals 
(relaxation pattern) are obtained according to the five-point first- 
and second-order differentiation coefficients in the left side of equa- 
tion (26). These coefficients for a 0.00001 change in the \fr-value at 
the grid point are given in table I. These and sim i lar coefficients 
for other value s of change in listed on cards were found to be con- 
venient in calculation. As a combination of checking and time saving, 
it is found convenient to relax a given set of residuals by using only 
the central three or five major coefficients at first and, when the 
residual is small enough, to compute for the total changes in i|r made 
at each grid point the correct resultant residual everywhere, and then 
to relax the residuals further with all coefficients. The technique of 
overrelaxing or underrelaxing and line-relaxing is very helpful . In 
most cases, the relaxation of 385 points (most relaxation is done in 
the blade region) is completed within 16 hours by hand computation. 

The greater portion of work for the present problem is obtaining the 
residuals, which takes about 40 and 60 hours, by hand computation, for 
the incompressible and the compressible solutions, respectively. In 
some of the later relaxation solutions obtained herein, this work of 
computation of residual is done on the IBM 604 Calculating Punch 
and an IBM Card Programmed Electronic Calculator (hereinafter called 
CPEC), except the computation of F r , which is more difficult to set 

up. In this way, the calculation of residuals for the compressible 
flow problem takes about 8 hours of machine time and a few hours of 
hand computation for F r . 
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Solution of Principal Equation by Matrix Method 

The 385 simultaneous algebraic equations (26) for the 385 interior 
grid points can be written in a compact matrix form by denoting the com- 
bined coefficients of fj's and V k ' s i n equation (26) at the point i 

by Ci and denoting the sum of N 1 and the product of the known boun- 
J . 

dary i|r -values and their corresponding coefficients by a 1 ! 


i|r^ = aA or 
3 


[c] {♦} - {a} 


(37) 


The matrix [c] is shown in figure 5 in terms of submatrices [E], [F], 
and [G]. Because there are seven interior grid points in the radial 
direction, these submatrices are of the order seven, and because of 
the cylindrical walls and the use of uniform spacing and the same dif- 
ferentiation formula, they repeat regularly along the diagonal of [c] . 

All the other submatrices are zero. The elements of [E], [F],[G] are 
given in figure 6 in terms of the grid spacings 5^ and 5*. In QfJ 
and [g] there are nonzero elements only on the diagonal. With the 
use of 5? and 5 * equal to 0.05 and 0.025, respectively, in order to 
get the elements of the matrices in short round numbers, the following 
multiplications are made: 


Row of [C] 

Value 
of r* 

Multipli- 

cation 

factor 

1, 8, 15, ... 379 

0.95 

0.57 

2, 9, 16, ... 380 

.90 

.54 

3, 10, 17, ... 381 

.85 

.51 

4, 11, 18, ... 382 

.80 

.48 

5, 12, 19, ... 383 

.75 

.45 

6, 13, 20, ... 384 

.70 

.42 

7, 14, 21, ... 385 

.65 

.39 


The resultant submatrices are shown in figure 7 . The same row multi- 
pliers are also applied to the elements of before solving for 

In the solution of' equation (37), because of the aforementioned 
special nature of [C], it is best to factor Qc] into two triangular 
matrices which also have not more than 15 nonzero elements in a row 
running alongside the diagonal (references 6, and 10 to 13). Thus 
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[C] = CL] [u] (38a) 

where [u] has elements which are unity along the diagonal, and equa- 
tion (37) becomes 

[L] [u] {+} = {a} (38b) 

The solution of {V} is then obtained by forward and backward substi- 
tution processes as follows: Let 

{3} = [u] {+} (39) 

Then 

[L] {p} = {a} (40) 

Solve | from equation (40) and then from equation (39). 

This matrix solution for the present investigation was made on an 
IBM CPEC and on a UNIVAC . Nine digits are used on the CPEC and eleven 
digits on the UNIVAC. At no place in the domain do the results 
obtained for {\|rj differ by more than 5 in the sixth digit. 

After the solution of is obtained in each cycle , it is sub- 

stituted into equation (37) for an over -all check. The residual found 
at any interior grid point is less than 1 in the eighth digit. Thus, 
the residual at every interior grid point as found in each cycle of 
calculation is reduced to practically zero with reference to the accu- 
racy of the grid size chosen. A comparison between the relaxation and 
the matrix methods showed that more accurate results were obtained 
with the matrix method in a smaller amount of total man- and machine - 
hours. 


RESULTS AND DISCUSSION 
Incompressible Solutions 

For adiabatic flow of a nonviscous gas with uniform s and H at 
the inlet, the deviation of gas flow from that on cylindrical surfaces 
in the incompressible case is due entirely to the term containing Fj* 
in equation (5). In order to study this effect, solutions are obtained 
for the two types of loading as given by equations (14) to (21). The 
one with uniform loading is designated case A, and the one with non- 
uniform loading, case B. 
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Convergence of solution. - Similar approaches are used in obtaining 
the solution for both cases. The calculation for case A is started with 
an assumption of straight cylindrical flow and successively corrected 
until the solution converges. The calculation for case B is started 
from an approximate solution of case A. The initial calculations in 
both cases are made by the relaxation method with 6 r = 0.1 and 

5* = 0.025. The values 5 r = 0.05 and 6 Z = 0.025 are used in later 
calculations. In each cycle, the \|r -values are difficult to improve 
beyond the fifth decimal place. 

The final relaxation solution is further improved by the matrix 
method. The \|r -values obtained in each matrix solution are resubstituted 
into equation (26) and the residuals calculated at any point are found to 
be less than 1 in the eighth decimal place. In the fourth matrix solution 
all values are sufficiently converged for the grid size and differentia- 
tion formula chosen. As an indication of convergence of the solution, 
the successive values of i|r at the mean radius r* = 0.8 in the matrix 
solutions of case B are shown in table II. The change in in the 

last cycle is less than 3 in the fifth significant figure or less than 
0.003 percent. The variations of i|r*, 3\|r*/ar*, 3\jr*/az*, F-J$r*, F^, F*, 
and N* at two points z* = 0.10 and 0.25 at the same radius 
obtained in the matrix solutions as well as the last three relaxation 
solutions are shown in table III. All results given in the following 
paragraphs are based on the last matrix solution and are the values on 
the mean stream surface as indicated in figure 4. 

Variation of F components and shape of streamline . - The calcula- 
tion is first started by using F r determined by the approx i mate for- 
mula (30) • It was later refined by using equations (7) to (9). The 
variations of F*r * and F* with z * at several radii are shown in 

figures 8(a) and 8(b), respectively. The variation of F*r* is very 
similar to that of s(V*r*)/8z* shown in figure 3^ being modified 
only by the variation in V*. 

The final values of p^. at the same radii are shown in fig- 
ure 8(c). The starting values determined by the approximate for- 
mulas (31) to (36) are also shown in the same figure and are seen to 
give reasonably good approximate values. The magnitude of F r i-s seen 
to be of the order of a quarter of F u (figs. 8(a) and 8(c)). This 
large value of F r is due to a combination of the large deflection of 

the gas in passing the turbine blading and the large radial twist of the 
stream surface of this type of velocity diagram and influences the flow 
distributions significantly. 
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Because of this F r -distribution, the nonhomogeneous term N in 
equation (5) or (26) is positive upstream of the radial element in the 
stator and downstream of the radial element in the rotor and is negative 
between them. The term N also increases in magnitude toward the hub. 
This distribution of N requires a general increase in \Jr at any point 
inside the blade row from its inlet value at the same radius and results 
in a streamline shape as shown in figure 9. In both cases A and B, the 
gas flows radially inward in the stator and outward in the rotor. The 
difference between the two cases is rather small. 

Variation of radial and axial velocities . - The variation of radial 
velocity with respect to z at several radii is shown in figure 10(a) . 

In both cases A and B, this velocity has a minimum and a maximum about 
in the middle of the stator and the rotor, respectively. It is practi- 
cally zero about 1 axial chord upstream and downstream. The radial loca- 
tion of the largest radial velocity occurs around r* = 0.75. 

The variation of axial velocity at five radii is shown in fig- 
ure 10(b). Its deviation from the inlet value is rather large along 
the casing and the hub , especially in the space between the two blade 
rows . 

/ 

Condition in plane normal to turbine axis and between stator and 
rotor . - In a plane normal to the turbine axis and between the stator 
and the rotor , the axial velocity is seen to decrease with an increase 
in the radius (fig. 10(h)). For radially uniform H and s, the equa- 
tion of motion in the radial direction is 


■3z 3r 


F 


r 


Because Fp is zero in the space between the stator and the rotor 
(fig. 8(c)), the radial velocity must decrease with an increase in z 
in order to satisfy this equation. This negative slope in V r with 
respect to z is clearly shown in case B by the values of radial 
velocity obtained at the regular grid points at z*= 0.15, 0.175, 
and 0.20 (fig* 10(a)). For case A, where V r is nearly equal to zero 
around z*= 0.15, this negative slope is apparent after the values 
of V r at z* = 0.1625 and 0.1875 are computed by using the three - 

point differentiation formula and the values of i|r at the regular grid 
points . 

It may be noted that this oscillatory variation of velocities in 
the space between the stator and the rotor is entirely due to the speci- 
fied variation of V u r in that space. Because V u r is specified to 
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be constant in that space (fig. 4), F u and F r on the stream surface 

are equal to zero in that space according to equations (7) and (9). 
These values are approached when the number of blades of the stator and 
rotor approaches infinity. In the actual case of a stream surface 
between two adjacent blades a finite distance apart, the streamlines 
turn for a short distance both upstream and downstream of the stator 
and the rotor (see, for example, reference 7), and consequently, both 
F r and the velocities would vary smoothly in the space between the 

stator and the rotor. 

Contours of constant velocity . - A contour plot of constant values 
of absolute velocity in the stator and relative velocity in the rotor 
for case B is given in figure 11. The maximum velocities leaving the 
stator and entering the rotor are, respectively, 1.77 and 1.25 times 
the rotor tip speed. 


Compressible Solution Neglecting F r 

A compressible solution for an inlet Mach number of 0.308 is first 
obtained by neglecting F r in the principal equation (case C). This 
solution is found in order to see the effect of the compressibility 
term in equation (5) alone and to see the error involved in neglecting 
F r for compressible flow. Nonuniform loading is used. 


The deviation of flow in this solution from cylindrical flow is 
due to the density term in N, which is mainly determined by the product 

fr fr ) ' ^his generally positive value of N inside the stator and 

the rotor requires a general decrease in \|r at any point inside the 
stator and the rotor from its inlet value at the same radius, thereby 
resulting in a streamline shape as shown in figure 12. 


In this solution, the starting values for the relaxation solution 
are obtained by using the values for the gas state at z* = 0.175 
obtained in reference 2 and by assuming that the streamline shape varies 
as a simple sine curve in the meridional plane. The solutions converged 
quickly. The successive values of \|r*, 8i|r*/3r*, 9 i|r*/ 3 z*, p*, and N* 
at r* = 0.8 obtained in the matrix and relaxation solutions are given 
in tables IV and V. The change in \|r* in the last cycle is less than 
2 in the fifth significant figure or about 0.002 percent. 

The variations of V*, V*, and p* at several radii obtained in 

the final matrix solution are shown for case C in figures 13(a), 13(b), 
and 14, respectively (the curves in those figures for cases D and E 
will be discussed in the following section) . Because of the decrease 
in density across the turbine, the axial velocity rises to a higher 
value downstream of the rotor (fig. 13(b)). 
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Compressible Solutions Including F r 

Solutions are next obtained for compressible flow by considering 
simultaneously the effects of compressibility and of F r . As in case C, 
nonuniform loading is used. Three axial positions of the radial ele- 
ment of the stream surface in the stator (z* = 0.0625, 0, and 0.15 for 

cases D, E, and F, respectively; see fig. 15) are considered in order 
to investigate the possibility of minimizing the effect of radial flow 
by the choice of this position. (Conceivably, the stream surface within 
the stator may have no radial element, its inclination with respect to 
a radial line being arbitrarily chosen to produce a certain desired 
effect . ) Because of the very large value of F r and the accompanying 

large radial flow, the solution for case F was stopped at an early 
stage , and only the value of Fr obtained in the early relaxation solu- 
tion and that obtained by the approximate formula (30) are given. The 
other two cases were carried further by the relaxation method and 
checked by two to three matrix cycles. Although the solutions are not 
so far converged as those in cases A, B, and C, they are accurate enough 
for ordinary purposes . The successive values of i|r and other pertinent 
variables at a number of typical points for case D are given in 
tables VI and VII in a manner similar to that for cases B and C. 

Nonlinearity of compressible solution . - In the calculation of 
case D, the values of the later relaxation solutions of case C were used 
as starting values. The results soon made clear, however, that quicker 
convergence would have been obtained if the calculation had been started 
from cylindrical flow. Because of the nonlinear nature of the principal 
equation, the complete solution for compressible flow for case D cannot 
be estimated from the solutions obtained in cases B and C. For the two 
contributions in the nonhomogeneous term N, the first term due to com- 
pressibility remains about the same as that in case C, but the second 
term, containing F r and p, is greatly changed because of the change 
in density. The distribution of the resultant values of N (and con- 
sequently the ^-distribution) is therefore quite different from that 
in case B or C or that obtained by directly combining the N-values in 
the two cases. 

Variation of F components . - The variations of F u r, F z , and F r 

are plotted as before in figures 16(a), 16(b), and 16(c), respectively. 
The variations of F u r and F z are similar to those in case B, which 
has the same type of loading. 

A comparison of the F r -distribution in the stator for the three 
values of z* chosen shows that alining the stream surface at the 
stator exit (z* = 0.15) causes very large negative values of F r in 
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the stator, which will produce large amounts of radial flow. . For 

z* = 0.0625, F r is negative for the first portion and positive for the 

second portion. For z$ = 0, F r is positive everywhere in the stator. 
The approximate formula (30) for Fr still compared very well with the 
final values. 

Shape of streamline. - The meridional projection of the streamlines 
obtained for cases D and E is shown in figure 17 . In case D, the gas 
flows radially outward and then inward in the stator, whereas in case E 
the gas flows radially inward and then outward in the s oator . The 
streamline in case D deviates from a cylindrical surface more than that 
in case E at the leading edge of the stator but less at the leading edge 
of the rotor. 

Variation of velocities and density . - The variations of V r and 
Vz in cases D and E at several radii are shown in figure 13. The radial 
velocities in these cases are of about the same order of magnitude as 
the velocity in case C but have more complicated shapes (fig. 13(a)). 
Figure 13(b) shows that there are significant differences between th$ 
values V z of the three cases C, D, and E, especially at the casing 
and the hub. These differences are due almost entirely to the different 
shapes in the streamlines because the differences in density among the 
three cases are very small (fig. 14) . 

Contour plots of static pressure and Mach number . - Because of the 
low value of the hub-tip radius ratio, there is^an over-all static - 
pressure rise across the rotor for values of r below 0.7 (fig. 18). 

An even greater static -pressure rise accompanied by a subsequent reexpan- 
sion to the exit pressure occurs within the rotor-blade rowj this 
increase in pressure rise is a result of the flow-area increase associ- 
ated with the assumption that B is equal to 1. 

The contours of constant Mach number (absolute Mach number in 
stator and relative Mach number in rotor) of cases D and E are shown in 
figure 19. The maximum Mach numbers at the stator exit are 0.982 and 
0.993 for cases D and E, respectively. If the effects of radial flow 
are completely ignored (simplified radial equilibrium), the maximum 
Mach number at the stator exit is only 0.927. 

Radial variation of axial velocity ahead of rotor . - In the 
girnpl ified— radial— equilibrium calculation, the axial velocity does not 
vary radially for free -vortex-type blading. The present solutions 
show, however, that there is considerable radial variation in axial 
velocity for both incompressible and compressible flows. An error in 
the axial velocity of the gas entering a blade row produces an error in 
Mach number and angle of attack, and thereby the range of efficient 
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operation is reduced. The radial variation of axial velocity ahead of 
the rotor (at z* = 0.175) is plotted for cases A to E in figure 20 in 
which the variation obtained in reference 2 is also plotted for com- 
parison. In all cases except simplified radial equilibrium, the axial 
velocity decreases with an increase in radius, and the rate of decrease 
is greatest at the hub. The effect of moving the radial element of the 
stream surface in the stator from the stator entrance (case E) to the 
midaxial position in the stator (case D) was to decrease the radial 
gradient in axial velocity at the rotor entrance. The approximate solu- 
tion of reference 2 is very close to the solution of case D. 

For incompressible flow, this variation is due to the shapes of 
the streamlines (fig. 9), as influenced by the Fr-distributions 
(fig. 8(c)), whereas for the compressible flow this variation is due 
mainly to the increasingly larger drop in density across the stator 
toward the hub. 


SUMMARY OF RESULTS 

A method recently developed for determining the steady flow of a 
nonviscous compressible fluid along a relative stream surface between 
adjacent blades in a turbomachine was applied to investigate the sub- 
sonic through flow in a single-stage axial-flow gas turbine. A free- 
vortex type of variation in tangential velocity was prescribed along 
the stream surface. Cylindrical bounding walls were specified in order 
to avoid radial flow at the walls. The flow variations on the stream 
surface for incompressible and compressible flows were obtained by using 
the relaxation method with hand computation and the matrix method on 
both an IBM Card Programmed Electronic Calculator and a UNIVAC. 

In all solutions considered in the present investigation, conver- 
gence was obtained without difficulty. A comparison between the relax- 
ation and the matrix methods showed that more accurate results were 
obtained with the matrix method in a shorter interval of time. The 
results of these highly accurate through-flow calculations form a basis 
for evaluation of simpler but less accurate methods. 

For incompressible flow, considerable radial flow was obtained 
because of the circumferential pressure gradient of the gas and the 
radial twist of the stream surface. The gas flowed radially inward in 
the stator and radially outward in the rotor. This radial flow resulted 
in a large negative radial gradient of axial velocity in the space 
between the stator- and rotor-blade rows. 

For compressible flow, the compressibility of the gas and the 
radial twist of the stream surface had equally important effects on 
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•the flow distribution and the nonlinear behavior of the principal equa- 
tion defining the flow was quite evident. The shape, or twist, of the 
stream surface and thus the amount of radial flow were sensitive to the 
axial location of the radial element of the stream surface in the stator. 
The largest radial flow occurred when this radial element was located at 
the exit of the stator. The effect of moving the radial element of the 
stream surface in the stator from the stator entrance to the midaxial 
position in the stator was to increase the deviation from cylindrical 
flow ahead of the stator and to reduce the deviation ahead of the rotor; 
the radial gradient in axial velocity at the rotor entrance was 
decreased. When the radial element of the stream surface in the rotor 
was near the midaxial position in the rotor, a large negative gradient 
of axial velocity was observed in all cases ahead of the rotor. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 18, 1952 
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TABLE I - RELAXATION COEFFICIENTS 


r 1 

r t 

For 0.00001 change of f at grid point (r^z^), 
change residual at (r,z) by 

r 

r t 

z 

z 1 - 25 z 

z 1 - 5 Z 

i 

z 

z 1 + 6 Z 

z 1 + 26 z 

0.95 

0.95 

.90 

.85 

-0.00133 

0.02133 

-0.04684 

.00519 

-.00031 

0.02133 

-0.00133 

0.90 

0.95 

.90 

.85 

.80 

-0.00133 

0.02133 

0.00232 

-.05000 

.00518 

-.00031 

0.02133 

-0.00133 

0.85 

0.95 

.90 

.85 

.80 

.75 

-0.00133 

0.02133 

0.00123 

.00548 

-.05000 

.00517 

-.00031 

0.02133 

-0.00133 

0.80 

0.95 

.90 

.85 

.80 

.75 

.70 

.65 

-0.00133 

0.02133 

-0.00032 

-.00035 

.00549 

-.06000 

.00516 

-.00031 

-.00036 

0.02133 

-0.00133 

0.75 

0.85 

.80 

.75 

.70 

.65 

-0.00133 

0.02133 

-0.00035 

.00550 

-.05000 

.00514 

.00149 

0.02133 

-0.00133 

0.70 

0.80 

.75 

.70 

.65 

-0.00133 

0.02133 

-0.00035 

.00551 

-.05000 

.00154 

0.02' 133 

-0.00133 

0.65 

0.75 

.70 

.65 

-0.00133 

0.02133 

-0.00036 

.00552 

-.04641 

0.02133 

-0.00133 
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TABLE II - SUCCESSIVE VALUES OF i|r IN MATRIX SOLUTIONS AT 

r* = 0.8 FOR CASE B 


* 

z 

Starting 

value 

Matrix solutions 

First 

Second 

Third 

Fourth 

CFEC 

UNIVAC 

-0.500 

0.08648 

0.086479 

0.086478 

0.086478 

0.086478 

0.086478 

-.400 

. 08648 

. 086475 

.086468 

.086468 

.086468 

.086468 

-.300 

.08648 

.086468 

.086450 

.086450 

.086449 

.086449 

-.200 

.08646 

.086448 

.086415 

.086415 

.086413 

.086412 

-.100 

.08639 

.086395 

.086354 

.086354 

.086350 

.086349 

-.050 

. 08636 

.086371 

.086330 

.086330 

.086324 

.086322 

-.025 

. 08638 

.086382 

.086333 

.086333 

.086327 

.086325 

0 

.08647 

.086437 

.086381 

.086381 

.086374 

.086372 

.025 

.08695 

.086879 

.086811 

.086811 

.086805 

.086803 

.050 

.08775 

.087673 

.087591 

.087591 

.087587 

.087585 

.075 

. 08882 

.088650 

.088550 

.088550 

.088550 

.088547 

.100 

.08976 

.089536 

.089415 

.089415 

.089419 

.089416 

.125 

. 09042 

.090145 

.089996 

.089996 

.090005 

.090002 

.150 

.09086 

.090522 

.090338 

.090338 

.090352 

.090351 

.175 

.09139 

.090972 

.090746 

.090746 

.090769 

.090769 

.200 

.09187 

.091458 

.091187 

.091187 

.091220 

.091221 

.225 

.09152 

.091119 

.090809 

.090810 

. 090834 

. 090835 

.250 

.09073 

.090336 

.089989 

.089989 

.090004 

.090007 

.275 

.08982 

.089464 

.089109 

.089107 

.089106 

.089108 

.300 

.08908 

.088756 

.088414 

.088414 

.088398 

.088400 

.325 

.08858 

.088284 

.087960 

.087960 

.087942 

.087944 

.350 

.08824 

.087981 

.087677 

.087677 

.087666 

.087668 

.375 

.08794 

.087728 

.087454 

.087454 

.087449 

.087450 

.400 

.08765 

.087515 

.087273 

.087273 

.087272 

.087273 

.450 

.08723 

.087185 

.087008 

.087008 

.087010 

.087011 

.550 

.08676 

.086800 

.086716 

.086716 

.086718 

.086718 

.650 

.08656 

.086622 

.086584 

.086584 

.086585 

.086585 

.750 

.08649 

.086537 

.086522 

.086522 

.086522 

.086522 

.850 

.08648 

.086489 

.086486 

.086486 

.086486 

.086486 


CM 

DO 


TABLE III - SUCCESSIVE VALUES OF \(r, F u r, F z , F r , AND N 

AT r* = 0.8 FOR CASE B 


z* 

Solution 

5* 

&z 

** 

37* 

3z * 

F* r* 

F* 

r z 

F* 

r r 

N * 

0.10 

Relaxation 

0.10 

0.025 

0.08987 

0.49248 

0.03500 

3.6851 

-7.2194 

0.49755 

0.58395 



.10 

.025 

.08982 

.49101 

.03563 

3.6741 

-7.2313 

.50430 

.59364 



.05 

.025 

.08976 

.49217 

.03238 

3.6828 

-7.2353 

.51134 

.60051 





a ( .08975) 








Matrix 

.05 

.025 

.0895364 

.49377 

.03050 

3.6948 

-7.2374 

.51534 

.60325 



.05 

.025 

.0894151 

.49266 

.02940 

3.6865 

-7.2386 

.51757 

.60723 



.05 

.025 

.0894189 

.49270 

.02959 

3.6867 

-7.2380 

.51589 

.60521 



.05 

.025 

.0894185 

.49269 

.02959 





0.25 

Relaxation 

0.10 

0.025 

0.09122 

0.48654 

-0.03647 

-5.4610 

-1.4682 

-0.19868 

0.23604 



.10 

.025 

.09098 

.49268 

-.03647 

-5.5230 

-1.4432 

-.20870 

.24484 



.05 

.025 

.09073 

.49292 

-.03603 

-5.5327 

-1.4453 

-.20328 

.23837 





a ( .09064) 








Matrix 

.05 

.025 

.0903360 

.49648 

-.03527 

-5.5726 

-1.4367 

-.20638 

.24837 



.05 

.025 

.0899893 

.49191 

-.03599 

-5.5213 

-1.4371 

-.20682 

.24302 

1 


.05 

.025 

.0900044 

.49235 

-.03667 

-5.5263 

-1.4376 

-.20674 

.24271 



.05 

.025 

.0900065 

.49243 

-.03665 






a Improved value By relaxation method and directly comparable to next value obtained by matrix 
method . 
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TABLE IV - SUCCESSIVE VALUES OF i|r IN MATRIX 
SOLUTIONS AT r* = 0.8 FOR CASE C 


* 

z 

Starting 

value 

Matrix solutions 

First 

Second 

Third 

CEEC 

UNIVAC 

-0.500 

0.08642 

0.086475 

0.086474 

0.086474 

0.086474 

-.400 

.08642 

.086451 

.086447 

.086447 

.086448 

-.300 

.08642 

.086399 

.086396 

.086396 

.086396 

-.200 

.08627 

.086285 

.086282 

.086283 

.086283 

-.100 

.08600 

.086024 

.086023 

.086023 

.086023 

-.050 

.08577 

.085789 

.085789 

.085789 

.085788 

-.025 

.08562 

.085631 

.085632 

.085632 

.085631 

0 

.08543 

.085440 

.085440 

.085441 

.085439 

.025 

.08520 

.085211 

.085210 

.085210 

.085210 

.050 

.08493 

.084949 

.084947 

.084947 

.084948 

.075 

.08467 

.084679 

.084676 

.084676 

.084676 

.100 

.08438 

.084439 

.084433 

.084433 

.084434 

.125 

.08417 

.084268 

.084258 

. 084258 

.084259 

.150 

.08414 

.084190 

.084179 

.084180 

.084181 

.175 

.08417 

.084218 

.084207 

.084207 

.084208 

.200 

: 08429 

.084352 

.084342 

.084342 

.084343 

.225 

.08454 

.084594 

.084585 

.084585 

.084586 

.250 

.08486 

.084884 

.084877 

.084877 

.084878 

.275 

.08515 

.085172 

.085166 

.085166 

.085167 

.300 

.08541 

.085424 

.085419 

.085419 

.085420 

.325 

.08562 

.085630 

.085626 

.085627 

.085627 

.350 

.08578 

.085796 

.085793 

.085793 

.085794 

.375 

.08591 

.085930 

.085928 

.085928 

.085928 

.400 

.08602 

.086040 

.086037 

.086037 

.086037 

.450 

.08618 

.086199 

.086195 

.086195 

.086195 

.550 

.08628 

.086370 

.086363 

.086363 

.086363 

.650 

.08644 

.086437 

.086432 

.086432 

.086433 

.750 

.08648 

.086463 

.086462 

.086462 

.086463 

.850 

.08648 

.086477 

.086477 

.086477 

.086477 


OJ 

0 ^ 


TABLE V - SUCCESSIVE VALUES OF 

r* = 0.8 



9\|r 

9z* 


FOR CASE C 


p, AND N AT 


* 

z 

Solution 

K 

5* 

z 


9ijf * 

9r* 

dZ* 

* 

P 

N* 

0.10 

Relaxation 

0.05 

0.025 

0.08442 

0.49085 

-0.01037 

0.82013 

0.19948 



.05 

.025 

.08438 

.49592 

-.01070 

.81871 

.20101 


Matrix 

.05 

.025 

.0844393 

.49552 

-.00840 

.81884 

. 19747 



.05 

. .025 

.0844328 

.49602 

-.00858 

.81870 

. 19748 



.05 

.025 

.0844335 

.49582 

-.00856 

.81876 


0.25 

Relaxation 

0.05 

0.025 

0.08486 

0.49087 

0.01253 

0.78938 

0.06684 



.05 

.025 

.08486 

.49638 

.01253 

.78771 

.06538 


Matrix 

.05 

.025 

.0848845 

.49512 

.01190 

.78810 

.06522 



.05 

.025 

.0848769 

.49535 

.01191 

.78802 

.06523 



.05 

.025 

.0848786 

.49528 

.01190 

.78804 

.06550 
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TABLE VI - SUCCESSIVE VALUES OF t IN MATRIX 
SOLUTIONS AT r* = 0.8 FOR CASE D 


z* 

Starting 

value 

Matrix solutions 

First 

Second 

Third 

CEEC 

UNIVAC 

-0.500 

0.08648 

0.086470 

0.086470 

0.086470 

0.086470 

-.400 

.08642 

.086422 

.086422 

.086419 

.086420 

-.300 

.08628 

.086330 

. 086330 

.086324 

.086324 

-.200 

.08597 

.086123 

.086123 

.086114 

.086114 

-.100 

.08539 

.085657 

.085657 

.085653 

.085651 

-.050 

.08497 

.085264 

.085264 

.085265 

.085262 

-.025 

.08474 

.085020 

.085021 

.085026 

.085023 

0 

.08452 

.084773 

.084773 

.084782 

.084778 

.025 

.08466 

.084860 

.084861 

.084872 

.084868 

.050 

.08511 

.085238 

.085239 

.085254 

.085249 

.075 

.08568 

.085730 

.085731 

.085753 

.085747 

.100 

.08614 

.086114 

.086114 

.086148 

.086140 

.125 

.08641 

.086312 

.086312 

.086358 

.086350 

.150 

.08660 

.086433 

.086433 

.086492 

.086483 

.175 

.08694 

.086713 

.086713 

.086785 

.086775 

.200 

.08741 

.087137 

.087137 

.087219 

.087207 

.225 

.08744 

.087124 

.087124 

.087220 

.087206 

.250 

.08720 

.086859 

.086859 

.086971 

.086953 

.275 

.08691 

.086558 

. 086558 

.086682 

.086661 

.300 

.08671 

.086361 

.086361 

.086498 

.086474 

.325 

.08662 

.086313 

.086313 

.086432 

.086407 

.350 

.08659 

.086331 

.086331 

.086433 

.086409 

.375 

.08656 

.086352 

.086353 

.086440 

.086417 

.400 

.08653 

.086373 

.086373 

.086448 

.086426 

.450 

.08650 

.086406 

.086406 

.086461 

.086443 

.550 

.08647 

.086447 

.086447 

.086474 

.086464 

.650 

.08648 

.086465 

.086465 

.086479 

.086473 

.750 

.08648 

.086474 

.086474 

.086480 

.086477 

.850 

.08648 

.086479 

.086479 

.086480 

.086479 


CM 

CD 


8 ^ 8 ^ 

TABLE VII - SUCCESSIVE VALUES OF § 7 , § 7 , P, AND N AT 

r* = 0.8 FOR CASE D 


z* 

Solution 

&r 



3V* 

3r* 

3^* 

3z* 

* 

P 

r r 

N* 

0.10 

Relaxation 

0.10 

0.025 

0.08616 

0.49474 

0.00888 

0.81905 

0.52876 

-0.28438 



.10 

.025 

.08609 

.49661 

.01593 

.81847 

.54071 

- .31213 



.05 

.025 

.08615 

.49675 

.01393 

.81845 

.54193 

-.30862 



.05 

.025 

.08614 

.49612 

.01450 

.81863 

.53724 

-.30484 





a (. 08614) 







Matrix 

.05 

.025 

.086114 

.49530 

.011534 

. 81889 

.54006 

-.30277 



.05 

.025 

.086148 

.49567 

.012020 

.81878 

.53905 

-.30187 



.05 

.025 

.086140 

.49564 

.011969 

.81879 



0.25 

Relaxation 

0.10 

0.025 

0.08737 

0.49501 

-0.00562 

0.78815 

0.20438 

-0.06956 



.10 

.025 

.08726 

.49521 

-.01177 

.78807 

.20388 

- .06909 



.05 

.025 

.08720 

.49493 

-.01180 

. 78815 

.20419 

-.06843 



.05 

.025 

.08720 

.49573 

-.01180 

.78791 

.20740 

-.06732 





a (. 08718) 







Matrix 

.05 

.025 

.086859 

.49450 

-.012527 

.78827 

.20320 

-.06944 



.05 

.025 

.086971 

.49483 

-.011951 

.78818 

.20331 

- .06983 



.05 

.025 

.086953 

.49488 

-.012092 

.78816 


1 


a Improved value by relaxation method and directly comparable to next value obtained by 
matrix method- 
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Figure 2. - Meridional section of gas turbine. 
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Figure 3. - Specified variation of tangential velocity V£r* 
with respect to axial distance z*. 
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Figure 4. - Grid system and relation between stream surface and reference meridional plane. 
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- Matrix [Cj expressed in terms of submatrices [E], [f], and [g]. 


G 

F 

E 


DO 


5 5 5 

1 ( 3 

1 1 

rH 

-f 

H | 

1 

0 

0 

0 

3a 2 0.95X6a 2b 2 

2a 2 0.95X2a 

3a 2 0.95X2a 

12a^ 0.95X12a 




4 2 

30 5 

4 ( 2 

1 1 

0 

0 

0 

3a 2 0 . 9X 3a 

12a 2 2b 2 

3a 2 0.9X3a 

12a 2 0.9x12a 




1 i 1 

4 2 

50 5 

4 . 2 

1 1 

0 

0 

12a 2 0.85X12a 

3a 2 0.85X3a 

12a 2 2b 2 

3a 2 0.85X3a 

12a 2 0.85X12a 



0 

1 t 1 

4 2 

50 5 

4 ( 2 

1 1 

0 

12a 2 0.8X12a 

3a 2 0.8X3a 

12a 2 2b 2 

3a 2 0.8X3a 

12a 2 0.8X12a 



0 

1 f 1 

4 2 

50 5 

4 2 

1 1 

0 

12a 2 0.75X12a 

3a 2 0.75X3a 

"l2a 2 2b 2 

3a 2 0.75X3a 

”l2a 2 0.75X12a 


0 


1 f 1 

4 2 

50 5 

4 ( 2 

0 

0 

12a 2 0.7X12a 

3a 2 0.7X3a 

12a 2 2b 2 

3a 2 0.7 X3a 




1 ; 1 

1 ( 1 

1 3 

5 : 5 5 

0 

0 

0 

12a 2 0.65X12a 

3a 2 0.65X2a 

2a 2 0.65X2a 

3a 2 0.65X6a 2b : 




(a) Submatrix [E ]. 
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Figure 6. 

- Submatrices [E], 

[F], and [G] expressed in terms of 

grid spacings 5* 

and b* z » a = 5*; 

b = 6». 
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( c ) Submatrix [G]. 
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Figure 6. - Concluded. Submatrices [e], [f], and [g] expressed in terms of grid 
spacings 5* and 8J. a = &*; b = 6*. 
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(a) Submatrix [E]. 
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(b) Submatrix [ F] . ( c ) Submatrix [g] • 

Figure 7. - Submatrices [E], [F], and [G] after row multiplication. 
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(a) Fjr*. 


Figure 8. - Variation of 


F 


in incompressible solutions 
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M Fj. 

Figure 8. - Continued. Variation of F in incompressible solutions. 
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(c) F*. 

Figure 8. - Concluded. Variation of F in incompressible solutions. 
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Meridional projection of streamlines of incompressible solutions 
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(a) Radial velocity V?. 

Figure 10. - Variation of velocities in incompressible solutions. 
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(b) Axial velocity V£. 

Figure 10. - Concluded. Variation of velocities in incompressible solutions. 
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Stator 


Rotor 



Figure 11. - Contours of constant absolute velocity In stator and relative velocity In rotor 

In case B. 
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z * 

Figure 12. - Meridional projection of streamlines for compressible solution neglecting F r . 
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Figure 13. - Variation of velocities in compressible solutions 
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Figure 14. - Variation of density In compressible solutions. 
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Figure 15. - Axial position of radial element of stream surfaces. 
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Figure 16 


Variation of F in compressible solutions 


N * 
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(b) F* for cases D and E. 


■igure 16. - Continued. Variation of F in compressible solutions. 
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Solution 


Approximate 
formula (30) 


Case 

D (z* = 0.0625 in stator) 
E (z* = 0 in stator) 

F (zq - 0.15 in stator) 



(c) Complete and approximate solutions for F * for cases D, E, and F. 
Figure 16. - Concluded. Variation of F in compressible solutions. 




Figure 17. - Meridional projection of streamlines of compressible solutions. 
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Figure 18. 
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(b) Case E. 

Figure 18. - Concluded. Contours of constant pressure ratio p/Pt,1* 


2402 



NACA TN 2750 


63 


r* 



Figure 19. 


(a) Case D. 

Contours of constant Mach number. 
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Figure 20. - Radial variation of axial velocity at z* equal to 

0.175. 
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